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In time-dependent density-functional theory, exchange and correlation (xc) beyond the adiabatic 
local density approximation can be described in terms of viscoelastic stresses in the electron liquid. 
In the time domain, this leads to a velocity-dependent xc vector potential with a memory containing 
short- and long-range components. The resulting time-dependent Kohn-Sham formalism describes 
the dynamics of electronic systems including decoherence and relaxation. For the example of col- 
lective charge-density oscillations in a quantum well, we illustrate the xc memory effects, clarify the 
dissipation mechanism, and extract intersubband relaxation rates for weak and strong excitations. 
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Time-dependent density-functional theory (TDDFT) 
Q has become a popular tool for describing the dynam- 
ics of many-electron systems. The exact time-dependent 
exchange-correlation (xc) potential v xc [n](r,t) contains 
information about the previous history of the system, in- 
cluding its initial state However, almost all present 
applications of TDDFT employ the adiabatic approxima- 
tion for v xc [n](r, i), ignoring all functional dependence on 
past time-dependent densities n(r, t'), t' < t. The sim- 
plest example is the adiabatic local-density approxima- 
tion (ALDA): 
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where e xc (n) is the xc energy density of a homogeneous 
electron gas of density n. The neglect of retardation in 
ALDA implies frequency-independent and real xc kernels 
in linear response 0]. This approach has been widely 
used for calculating molecular excitation energies 0, . 

There have been several attempts to go beyond the 
ALDA H ©, 0, H H OH- Vignale and Kohn (VK) 
showed that a non-adiabatic local approximation re- 
quires the time-dependent current j(r, i) as basic vari- 
able, rather than the density n(r, t). This formalism was 
later cast in a physically more transparent form using 
the language of hydrodynamics [7], |8|, M, LL0( , where non- 
adiabatic xc effects appear as viscoelastic stresses in the 
electron liquid. 

To date, the VK formalism has been applied exclu- 
sively in frequency-dependent linear response. The first 
application was to calculate linewidths of intersubband 
plasmons in semiconductor quantum wells In the 

absence of disorder and phonon scattering, ALDA gives 
infinitely sharp plasmon lines. The VK formalism in- 
cludes damping due to electronic many-body effects, in 
good agreement with experimental linewidths [ill ]. Van 
Faassen et al. ^1 calculated static axial polarizabilities 
in molecular chains, which are greatly overestimated with 
ALDA. For many systems, VK led to an excellent agree- 
ment with ab initio quantum chemical results. Other 
recent studies applied the VK theory to calculate atomic 



and molecular excitation energies 0, Q . 

Beyond linear response, a wealth of interesting electron 
dynamics can be explored using time-dependent Kohn- 
Sham (TDKS) theory 0]. This letter presents an analy- 
sis and applications of the VK formalism in the time do- 
main. A striking consequence of the memory and velocity 
dependence of the VK xc potential is that it introduces 
retardation into TDKS theory, which in turn leads to de- 
coherence and energy relaxation. This will be illustrated 
for charge-density oscillations in quantum wells. 

Several time-dependent Schrodinger equations with 
dissipation have been proposed in the literature 0, Q, 
firL Il8| . using quantized classical frictional forces or other 
phenomenological assumptions. By contrast, the VK xc 
potential has a microscopic origin, and satisfies exact con- 
ditions such as the harmonic potential theorem 0, S 0| • 

In the presence of external scalar and vector potentials, 
v(r, t) and a(r, t), the TDKS equation is 
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where «h is the Hartree potential and a xc (r, t) is the xc 
vector potential. A non-adiabatic, nonlinear xc vector 
potential has been given in Ref. to within second 
order in the spatial derivatives as: 

1 a xc ,,(r, t) = -V iWx A c LDA (r, t) + £ Yl^iiM , (3 ) 
c n(r,t) 

j 

where the viscoelastic stress tensor er xc is defined in terms 
of the velocity field u(r, t) = j(r,i)/n(r, t): 
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The viscosity coefficients in Eq. (0J are defined as 
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and similar for £, where 
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f£ c and f£ c are the longitudinal and transverse frequency- 
dependent xc kernels of a homogeneous electron gas of 
density n 0, 0|. The apparent ambiguity in Eq. (JSJ) 
whether the density should be evaluated at t or t' is re- 
solved by noting that the difference involves higher gra- 
dient corrections. The same argument applies to the dif- 
ference between the instantaneous position r of a fluid 
element and its "retarded position" R 

In the following, we consider quantum well systems 
where all spatial dependence is along the z direction 
only. One can then transform the xc vector potential, Eq. 
©, into a scalar one: v xc (z,i) = u xc LDA (z,i) + v^.(z,t) 
(ALDA+M), with the memory part given by 
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Assuming that the system has been in the ground state 
(with zero velocity field) for t < 0, the zz component of 
the xc stress tensor becomes 



a xc , zz (z',t)= / Y(n(z',t),t-t')V z .u' z (z',t')dt' , (9) 
Jo 

where the memory kernel Y is given by 



Y(n, t-t') = - 7j(n, t-t') + C(n, t-t') . (10) 
3 

With the help of the Kramers-Kronig relations for f~ c we 
can express the memory kernel as follows: 
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with the xc shear modulus of the electron liquid |l9j 
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Fig. n shows the memory kernel Y for r s = 3, eval- 
uated with the Gross-Kohn (GK) and Qian-Vignale 
(QV) [l9j parametrizations for f^ c {u>). As shown in the 
inset, Y GK decreases exponentially over time. The falloff 
is very rapid for the highest densities, and the memory 
becomes more and more long-ranged for lower densities. 
Comparing the GK and QV parametrizations, one finds 
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FIG. 1: Memory kernel Y(n . t- t') for r s = 3, using the QV 
and GK parametrizations |3lll9|| for f^ c {uj). Inset: y GK for r s 
between 1 and 15, indicating exponential memory loss, with 
a longer-ranged memory for lower densities. 



a similar overall behavior for different values of r s , ex- 
cept that yQ v does not decrease monotonically with time 
but passes through a negative minimum, and then ap- 
proaches the finite limit 2/it xc /3 (/i xc — * for large r s 
[19j, while GK assume /Lt xc = throughout). 

To illustrate and analyze the xc memory effects beyond 
the ALDA, it is convenient to use a simple analytic model 
density in Eq. J5J) to evaluate v^.(z,t). The function 
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mimics the non-interacting density of a hard-wall quan- 
tum well, driven by an AC field of frequency lj, where 
\A\ < 1 to ensure that n(z, t) > 0. We take N s = 1 a^ 2 
and L = 10 ao such that r s ~ 1 in the center, and we sim- 
ulate a weakly driven case with to = 1 a.u. and A — 0.01. 

Fig. [2]shows a stroboscopic plot of v^.(z, t) during the 
4th cycle after switch-on of the time dependent model 
density. Four snapshots during the cycle are highlighted: 
density passes through equilibrium, turns around at the 
right wall, sloshes back, and hits the left wall. The re- 
action of v^. to these periodic density fluctuations is re- 
markable: it opposes the instantaneous current flow by 
periodically building up an S-like potential barrier in the 
central part of the well, trying to slow down the sloshing 
motion of the density. The large-amplitude fluctuations 
of close to the edges, on the other hand, have little 
overall impact since they occur in low-density regions. 

Both v^,' and u^'^ v have similar magnitude and 
shape, and lag behind the fluctuations of u A c LDA , giv- 
ing rise to a retardation of the ALDA+M xc potential. 
From the data in Fig. [2J v x ^' C ^ v has a quarter-cycle (tt/2) 
phase lag, which is the most efficient for retardation: the 
opposing potential barriers are largest at the instants of 
maximal current flow, and become flat whenever the den- 
sity hits the wall and turns. The good synchronization of 
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FIG. 2: Memory part of the xc potential [Eq. JSJ], evaluated 
for n(z, t) of Eq. 1131 in GK and QV parametrization, shown 
in a stroboscopic plot during the 4th cycle after switch-on. 
The heavy lines indicate equidistant snapshots. Compared to 
the ALDA fluctuations (bottom panel), i>^ ,C3V has a phase 
lag 7r/2, and f^' GK a little over 7i"/2. The density oscillation 
in the inset is drawn with enhanced amplitude for clarity. 



QV points to a delicate balance between the short- and 
long-range parts in the memory kernel Y"Q V (see Fig. Q). 
By comparison, the timing of v^.' GK is slightly off with a 
phase lag greater than tt/2. As a consequence, QV causes 
stronger damping than GK (see below), even though its 
potential barriers are lower. 

We now turn to a more realistic case and solve 
the TDKS equation for a 40-nm modulation-n-doped 
GaAs/Alo.3Gao.7As quantum well 0, with effec- 
tive mass m* = 0.067m, effective charge e* = e/VT3, 
well depth 257.6 meV, and electronic sheet density N s = 
1 x 10 11 cm~ 2 . The initial condition, obtained from the 
static Kohn-Sham equation in LDA, is the electronic 
ground state in the presence of a uniform electric field 
£ ("tilted" quantum well). At t = 0, the electric field 
is abruptly switched off, which leaves the quantum well 



0.0015 




-0.06 1 1 ' ' ' 1 

20 40 60 80 100 

t (a.u.) 

FIG. 3: Dipole moment d(t) of a 40-nm GaAs/Al .3Ga .7As 
quantum well with electron density 1 x 10 11 cm -2 , initially 
in a uniform electric held S\ = 0.01 mV/nm (top) and £2 = 
0.5mV/nm (bottom), which is abruptly switched off at t — 0. 
Dashed lines: ALDA. Full lines: ALDA+M (using QV). 

electrons in an excited state and triggers a collective 
charge-density oscillation. In the linear regime, this oscil- 
lation represents the so-called intersubband (ISB) plas- 
mon, but in the following we will also explore the non- 
linear, strongly excited regime. Fig. [3] shows the dipole 
moment d(t) = J z n(z,t)dz versus time (1 a.u. = 61 fs 
in GaAs) for two initial electric fields, £\ — 0.01 and 
£2 = 0.5 mV/nm, comparing ALDA with ALDA+M (us- 
ing QV; GK gives qualitatively similar results). 

Linear-response theory yields ISB plasmon fre- 
quencies walda = 0.9307, wqk = 0.9367, and wqv = 
0.9348 a.u., and line(half)widths r G K = 0.0061 and 
F Q v = 0.0078 a.u. (1 a.u. = 10.79 meV in GaAs). The 
weak-field case, £1, is perfectly reproduced by d{t) = 
do cos(walda0 in ALDA and d(t) = do cos(uiQvt)e~ rQvt 
in ALDA+M. This shows that v^. has dissipative and re- 
active components, giving an exponential decay of the 
dipole amplitude with characteristic time Tgy, and a 
small blueshift of the ISB plasmon frequency. 

The stronger field, £2, causes much richer dynamics. 
By projecting the initial wavefunction on field-free Kohn- 
Sham states, we find 10% and 0.1% initial occupation 
probabilities of the second and third subband levels (see 
Fig. 0] inset). Consequently, a more complex pattern 
emerges in the dipole oscillations. A spectral analysis 
shows that d(t) contains higher harmonics of the lowest 
(1-2) ISB plasmon, as well as sidebands resulting from 
a nonlinear coupling with the dipole-forbidden 1-3 ISB 
excitation (the crosstalk between the 1-2 and 1-3 modes 
is mediated through modulations of the TDKS effective 
potential). Due to their larger velocities, the higher- 
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FIG. 4: Dissipation of excitation energy E(t), for the quan- 
tum well of Fig. [jj] and initial fields £ between 0.01 and 1 
mV/nm as indicated, calculated with ALDA+M (using QV). 
Inset: initial occupation probabilities of the first three sub- 
bands. 



bution as a superposition of local plasmon modes, each 
subject to decay into multiple particle-hole excitations. 

In conclusion, we have given an explicit demonstration 
how memory effects introduce the element of intrinsic de- 
coherence and energy relaxation into TDKS theory. This 
represents an alternative viewpoint to the density-matrix 
approach |2flj . which, in its simplest form, describes dis- 
sipation through phenomenological decoherence and re- 
laxation times (known as T\ and T2 for 2-level systems) . 
A combination of the two approaches suggests itself as 
a powerful TDDFT tool to describe nonlinear electron 
dynamics in the presence of intrinsic and extrinsic (im- 
purities and disorder) dissipation mechanisms. 
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frequency spectral components in d(t) are more rapidly 
damped than the 1-2 ISB plasmon, as seen in Fig. 

The decoherence of the dipole oscillations is accom- 
panied by energy dissipation. We define the excitation 
energy per unit area, E(t), as the expectation value of 
the ALDA Hamiltonian with the ALDA+M wave func- 
tion, minus the ground-state energy. Fig. 0] shows E(t) 
scaled by the square of the initial electric field. Following 
the quadratic Stark effect, E{t)/£ 2 is independent of £ 
for small fields < O.OlmV/nm. For larger £, higher-order 
deviations from the quadratic Stark effect emerge. 

For small £ , the excitation energy decreases as E(t) = 
E(0)e~ 2rQvt (small steps are superimposed since the in- 
stantaneous dissipation rate depends on the oscillating 
velocity field of the ISB plasmon). For larger £ , one 
notes deviations from this simple behavior in the form of 
a more rapid initial decay and the appearance of a larger 
step structure. The origin for these steps is the nonlin- 
ear coupling between the 1-2 and 1-3 ISB plasmon modes 
discussed above, which generates sidebands around Wq 2 v 
whose frequency spacing ujq V — 2wq 2 v increases with £ 
(0.27, 0.39, and 0.52 a.u. for 0.1, 0.5, and 1 mV/nm). As 
long as these steps are not too pronounced, E(t) is well 
described by a biexponential model, with an additional 
fast channel accounting for the relaxation from higher 
subbands. The associated relaxation rate varies between 
0.017 and 0.021 a.u. for £ between 0.1 and 0.3 mV/nm, 
which is more than twice as fast as Fqv- 

Finally, we comment on the physical mechanism for 
energy dissipation. In the linear regime, the VK the- 
ory locally assumes a homogeneous electron gas subject 
to small periodic modulations. The frequency-dependent 
xc kernels then cause the decay of collective modes into 
multiple particle-hole excitations, even if Landau damp- 
ing is forbidden 0, 0, 13 • In the time domain, one can 
view the dynamics of an inhomogeneous electron distri- 
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